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We present the reconstruction of the Wigner function of a classical phase- 
sensitive state, a pulsed coherent state, by measurements of the distributions 
of detected-photons of the state displaced by a coherent probe field. By 
'. using a hybrid photodetector operated above its photon-resolving regime, 

. we obtained both the statistics at different values of the probe field and the 
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H ' 1. Introduction 

Quantum optical states are a valuable resource in view of the development of quantum pro- 
tocols as the many degrees of freedom of the radiation field can be exploited to enhance 
transmission, manipulation and storage of quantum information. To this aim, during the 
last decades many efforts have been devoted to the generation of suitable quantum states 
of light: besides the production of different kinds of squeezed states [1,2], conditional states 
endowed with non-Gaussian nature have been produced starting from nonlinear paramet- 
ric processes [3-7]. Among them phase-sensitive quantum states can be considered a good 
resource as they also carry phase information [8]. Obviously, to exploit such states reliable 
characterization strategies are required. The typical technique employed so far is quantum 
homodyne tomography (OHT) [9] operating either in spectral or in time domain [10]. OHT 



1 



was proposed in 1989 by Vogel and Risken [11] and is based on balanced homodyne detec- 
tion [12,13], where the interference of a signal field and a local oscillator (LO), a coherent 
field with variable phase, provides the values of the quadrature amplitude distributions of 
the signal as a function of the LO phase. In general the variation of the phase is obtained 
by means of piezoelectric transducers controlled by a feedback system. In order to overcome 
the unwanted fluctuations arising from the instabilities of the movable elements, quite re- 
cently high-frequency time-domain OHT has been proposed and implemented [14-16]. This 
technique can operate at the high repetition rate available from mode-locked lasers. Unfortu- 
nately, to obtain quantum states containing sizeable numbers of photons, low rep-rate laser 
sources must be used, and the acquisition of OHT data may become unpractically long and 
affected by uncontrollable instabilities. 

To overcome this problem, here we propose a self-consistent Wigner function reconstruction 
based on the interferometric scheme depicted in Fig. [3] [17]. The method consists in detect- 
ing the light exiting a beam-splitter that mixes the signal field to be characterized with 
a coherent probe field, whose amplitude, \a\, and phase, <f>, can be varied with continuity. 
We have already applied this method to the measurement of the Wigner function of some 
phase-insensitive fields [18] produced by a low rep-rate source (15 kHz); now we present a 
measurement on a phase- sensitive coherent field from the same source. At such a rate, the 
measurements took a long time and the problem of determining the phase in the presence of 
potential instabilities was solved by exploiting the linearity of the detector in a completely 
self-consistent analysis procedure. 

The paper is organized as follows: in Section [2] we present the theoretical background of 
the method, in Section [3] we present in detail experimental procedure and results and finally 
we conclude in Section HJ 

2. Theory 

The Wigner function of a single mode state can be written as [17] 



where p n (a) is the photon- number distribution of the field displaced by a, which is not di- 
rectly accessible by using real detectors having quantum efficiency r] < 1. 
However, if we are able to characterize the detector, knowing the statistics of detected pho- 
tons can be enough for a full characterization of the signal field. For instance, if we assume 
that the primary detection process is described by a Bernoullian convolution, we have that 
the detected-photon distribution, p m (/3), is given by [19] 




n=0 
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where (3 = Jr\a. is the detected amplitude of the displacement field. 

In principle, p m (P) ^ p n ( a ), that is the distribution of the photoelectrons is not necessarily 
the same as that of the photons. Nevertheless it can be demonstrated [20,21] that when the 
state can be obtained by combining non-squeezed gaussian fields, the statistics of detected 
photons remains the same. 

If Eq. (J2J) holds, we can write for the Wigner function of the detected photons an expression 
analogous to that in Eq. [18] 

c\ OO 

W(f3) = -Y.(- i rPm(f3) , (3) 
n m =o 

where 

WW) = ^(1^) / d 2 P'e-^-^W(P'/^j) . (4) 

is the Wigner function in the presence of losses [22] and W(f3/ \frj) = W(a) is that of the 
photons. Obviously, if p m (/3) has the same functional form as p n (a), also the Wigner function 
remains the same. 

As an axample, in the case of a coherent signal, the Wigner function for the detected state 
is 

W((3) = -ex P (-2\(3-(3 \ 2 ) , (5) 

where (3q is the complex amplitude of the considered coherent signal field. 

The statistics of detected photons can be obtained with some classes of detectors, including 
photomultipliers [23,24], cryogenic photon counters [25,26], multi-pixel photon counters [27], 
time-multiplexed detectors [28,29], multichannel fiber loop detectors [30]. Obtaining or not 
the statistics of detected photons depend on the photon-counting capability of the detector, 
a fact that limits the range of intensities that can be investigated. For all these reasons, the 
only experimental Wigner functions directly measured by photon counting are those reported 
in Ref. [31] for very low intensities. 

Recently we introduced a new method to analyze the output signal of linear detectors, 
that allowed us to reconstruct the detected photon distribution for a number of optical 
states without any a priori knowledge [20,32,33]. In our technique we consider a linear 
photodetector endowed with high gain and model the detection process through a Bernoullian 
convolution [34] and the overall amplification and conversion process through a very precise 
factor, which is taken as constant, 7 [33]. With these assumptions, we can link the statistics 
of the detector output (voltages in the present case) to that of the detected photons. In 
particular, we can write the mean value as v = 7m = 7777!, 77 being the overall detection 
efficiency, and the variance as o~ 2 (v) = 7V 2 (m) = r y 2 (a 2 (n) + r](l — ri)n). We can thus measure 
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the Fano factor of the output voltages, F v = a 2 (v)/v, at different values of 77. We observe 
that F v can be written in the very general form 



F v = —v + 7 
n 



(6) 



where all the dependence on the field under investigation is in the angular coefficient, Q/n, 
Q being the Mandel parameter. We can thus obtain 7 from a fit of experimental data. Once 
7 is evaluated, it is possible to find the photoelectron distribution by dividing the v output 
values by the experimental value of 7 and re-binning the data in unitary bins. Our method 
has the advantage of being self-consistent as the value of 7 is obtained from measurements 
on the same field under investigation [32] and does not require independent calibration. 

3. Experiment 

Aim of the experiment is to reconstruct the Wigner function of a phase-sensitive state (sig- 
nal), a coherent state, by using Eq. (JHJ). To reach this goal, the probability distribution 
of detected photons must be measured at different values of the complex displacement, 
a = \a\ exp(z0) (probe), to obtain p m (|/3| exp(i0) = y^|o:| exp(z0)). The experimental setup 
is sketched in Fig. [3j 

The light source was a frequency-doubled Q-switched Nd:YAG laser operating at 15 kHz rep- 



Nd:YAG 











spatial filter 



SGI 



BS 



ADC+PC 



BS k 

m 

probe p/ 



fx 



piezo-driver 



Fig. 1. (Color online) Experimental setup. HPD, hybrid photo-detector; BS, 
beam splitter; F, neutral density filter; Pz, piezoelectric movement; P, polar- 
izer; SGI, synchronous gated integrator. 



etition rate (Quanta System) and providing linearly polarized pulses at 532 nm of ~ 200 ns 
duration. The beam was spatially filtered and split into two parts serving as signal and 
probe fields. Intensity and phase of the probe field were adjusted by means of a variable 
neutral filter (F) and of a piezoelectric movement (Pz), respectively. Signal and probe were 
then mixed at a cube beam-splitter (BS in Fig. [3]). The field exiting one of the ports of the 
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BS was passed through a diverging lens and a portion of it was collected by a multimode 
optical fiber (600 /im core-diameter) and delivered to a hybrid photodiode module (HPD, 
H8236-40 with maximum quantum efficiency ~ 0.4 at 550 nm, Hamamatsu) operated above 
its photon-resolving regime. In fact, this detector is endowed with a limited photon-number 
resolving capability but is linear over a wide range of intensities. 

We sampled the phase space by changing amplitude and phase of the probe independently. 
In particular, we considered 46 values of probe intensity, |/3| 2 , set by adjusting the variable 
filter F, and 48 values of phase, 0, at each set by moving the Pz in steps. Before starting 
the piezoelectric movement, we made independent measurements of both signal and probe, 
and evaluated their mean values \(3q\ 2 and \/3\ 2 , respectively. At each one of the 46x48-values 
of the probe, we performed the calibration procedure described in Section El to obtain the 
value of 7. Each calibration run included 25 ?7-values, set by means of the polarizer P in 
Fig. [3] from the maximum value t] max = 0.31 down to almost zero. For each 77 we measured 
the output voltages, v, at 30000 single-shots, we evaluated the Fano factor, F v , and from its 
fit as a function of v we obtained the value of 7 according to Eq. (jBJ). By using this value 
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Fig. 2. (Color online) Columns: measurements taken at the indicated values 
of probe, |/3| 2 , and signal, |/3o| 2 , intensities. First row: Wigner function recon- 
structed according to Eq. as a function of the piezo step, i. Second row: 
mean value of the output voltages at r] max . Third row: experimental cosine of 
the interference pattern. Fourth row: phase values obtained from the cosine. 

of 7, we calculated the probability distributions p m (0) corresponding to r\ max in each run 
and used them to calculate W{jS) according to Eq. ([3]). The resulting values of the Wigner 
functions are plotted as a function of the step of the piezo, i, in the first row of Fig. [3j 
We observe that, due to unavoidable instabilities, the reproducibility of the movement of the 
piezo was not guaranteed. As a consequence, the Wigner functions thus reconstructed cannot 
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be directly interpreted. Moreover, although a large part of the experimental procedure was 
automated, the complete acquisition required very long measurement sessions, performed in 
different days. For all these reasons, we had to devise a strategy to recover the phase values 
directly from the measurement by taking advantage of the linearity of the detector. 
The procedure is outlined in the following. At each step of the piezo movement, we calcu- 
lated the mean value, v, of the most intense acquisition, corresponding to r] max (see second 
row in Fig. [3]). For regular steps of the piezo, we would expect a sinusoidal trend repro- 
ducing the interference pattern between signal and probe, but, in the absence of a feedback 
control on the piezo, its positioning was extremely irregular. Nevertheless, the pattern still 
represents an interference and thus we used the data to infer fringe visibility, and hence 
overlap of signal and probe, and phase difference between signal and probe. We calculated 
the interference visibility as V = (max(u) — min (?;))/ (max (y) + min(ir)) and the overlap as 
£ = V max /{2 — V max ) [35], where V max was obtained when the intensity of the probe was 
equal to that of the signal. In our system we had a limited visibility, V max = 0.87, and hence 
£ ~ 0.78. The visibility values for the data reported in Fig. [3] were V = 0.73, V = 0.77 and 
V = 0.86. To obtain the phase, 0, we evaluated a quantity that can be interpreted as a cosine 
of the phase, namely cos(0) = (v — (max(u) + min(w))/2)/(max(u) + min(u))/2 (see third 
row in Fig. [3]). By inverting the cosine we obtained the values of <fi displayed in the fourth 
row of the same figure. This method allows assigning the actual value of phase at each piezo 
position independent of the regularity and reproducibility of the movement. 
In Fig. [3] we present as dots the values of the experimental Wigner function, same data as 
in Fig. plotted as a function of the obtained <fi values. Together with the experimental 
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Fig. 3. (Color online) Dots: experimental reconstruction of sections of the 
Wigner function at fixed values of probe, |/9| 2 , and signal, |/^o| 2 , intensities and 
variable 0; dash-dotted lines: expected theoretical curves; full lines: theoretical 
curves corrected for the overlap of signal and probe. 

data, in Fig. [3] we also plot the expected theoretical curve evaluated from Eq. (JSJ) at the 
measured values \f3\ and \f3 Q \ (dash-dotted line). The apparent discrepancy between theory 
and experiment can be ascribed to the non-perfect overlap, £, of signal and probe: only a 
part of the measured state comes from the superposition while the remaining is simply a 
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spurious coherent contribution coming from a part of the probe. As a result, the measured 
distributions become the convolution of those expected and the poissonian distribution of 
the residual coherent probe [18]. In the present case, in which all the distributions are pois- 
sonian, no variation is visible in the p m (/3), while the theoretical Wigner function must be 
recalculated as the product of the Wigner function expected for a signal reduced by the 
overlap with the gaussian Wigner function of the residual field [31]: 

W(/3) = ^exp(-2(l-OI/5| 2 )exp(-2| v /e/3-/5 | 2 ) • (7) 

The theoretical Wigner functions, calculated from Eq. ([7]) with £ = 0.78 and \/3q\ = 0.70, 
\Po\ = 0.63 and |/3 | = 0.67, respectively, superimpose to the experimental data (see full lines 
in Fig. [3]). 

Fig. 4. (Color online) Main: experimental reconstruction of the Wigner function (dots) and 




theoretical prediction. Inset: contour plot of the experimental data. Finally, Fig. [3] shows the 
3D reconstruction of the Wigner function obtained by combining all the measurements at 
different (3 (dots). The theoretical surface was obtained by correcting the expected Wigner 
function in Eq. (jSJ), taking into account the overlap and setting the mean value of the signal 
intensity at \(3q\ = 0.67. For a quantitative estimation of the quality of the reconstruction, 
we evaluate the mean error e 

= Ek=iW {Pk) - W c * P (p k )]/N. The presented data yielded 
e = — 1.7 x 10 -3 . As we can see, the theoretical surface carefully fits the experimental data 
in spite of all the experimental instabilities, thus testifying the robustness of both setup and 
analysis method. 

4. Conclusions 

In conclusion, we implemented the reconstruction of the Wigner function of a pulsed, phase- 
sensitive, coherent field state containing a low, but non- negligible, number of photons (n = 
m/Vmax = 2.12). The technique takes advantage of the features of a hybrid photodetector, 
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which is both able to distinguish a single detected photon and endowed with linear response 
over a wide intensity range. The heart of our experiment is the self-consistent procedure 
employed to analyze the detector output that allowed us to establish the statistics of detected 
fields and to assign the values to the phase of the probe field. The technique results to be 
reliable in spite of the low rep-rate of the laser source that imposed long experimental sessions 
to collect the data. This opens the way to the characterization of non-classical states, such 
as conditional states generated starting from twin-beam and squeezed states containing a 
sizeable number of photons, that can only be generated by low-rep-rate sources. Moreover, 
the self-consistent calibration procedure, which in the present work was repeated at each of 
the 46 x 48 values of the probe, as expected yielded almost constant values of 7 once we had 
set the parameters of the acquisition chain. Such a stability in 7 should allow us to perform 
the calibration procedure only at the beginning of the measurement, thus shortening the 
overall measurement time. 
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Fig. 1. (Color online) Experimental setup. HPD, hybrid photo-detector; BS, beam splitter; 
F, neutral density filter; Pz, piezoelectric movement; P, polarizer; SGI, synchronous gated 
integrator. 
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Fig. 2. (Color online) Columns: measurements taken at the indicated values of probe, |/3| 2 , 
and signal, |/?o| 2 , intensities. First row: Wigner function reconstructed according to Eq. ([3]) as 
a function of the piezo step, i. Second row: mean value of the output voltages at r] max . Third 
row: experimental cosine of the interference pattern. Fourth row: phase values obtained from 
the cosine. 
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Fig. 3 (Color online) Dots: experimental reconstruction of sections of the Wigner function 
at fixed values of probe, |/3| 2 , and signal, \/3q\ 2 , intensities and variable 0; dash-dotted lines: 
expected theoretical curves; full lines: theoretical curves corrected for the overlap of signal 
and probe. 
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Fig. 4. (Color online) Main: experimental reconstruction of the Wigner function (dots) and 
theoretical prediction. Inset: contour plot of the experimental data. 
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